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Abstract 

I study some classes of RG flows in three dimensions that are classically conformal and have 
manifest g —* 1/g dualities. The RG flow interpolates between known (four-fermion, Wilson- 
Fischer, tp®) and new interacting fixed points. These models have two remarkable properties: 
i) the RG flow can be integrated for arbitrarily large values of the couplings g at each order of 
the 1/N expansion; ii) the duality symmetries are exact at each order of the 1/N expansion. I 
integrate the RG flow explicitly to the order 0(1/N), write correlators at the leading-log level 
and study the interpolation between the fixed points. I examine how duality is implemented in 
the regularized theory and verified in the results of this paper. 
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1 Introduction 



In three-dimensional quantum field theory it is often possible to derive non-perturbative results 
using a large N expansion. Typically, N is the number of scalar fields or fermions. To a 
given order of the 1/N expansion infinitely many graphs can be collected by means of simple 
resummations and the resulting number of effective graphs is finite. In some cases it is not 
necessary to assume that the couplings g of the theory are small and expand perturbatively in 
g, but it is possible to calculate the exact ^-dependence at each order of the 1/N expansion. 
Then, the computation of a finite number of effective graphs allows us to integrate the RG flow 
and study the interpolation between the fixed points. If the theories satisfy a duality g <-> 1/g, 
then this symmetry is implemented exactly at each order of the 1/N expansion. The purpose 
of this paper, which is an elaboration of the ideas of ref. [0, is to study a variety of models 
with these features. 

Specifically, I study two classes of flows. The first class of models involves K + 1 sets of Nj, 
I = 0, 1, ... if, fermions interacting by means of a cr-field, that is to say a dynamical source 
for a composite operator. The source a is called dynamical because the functional integral is 
performed also over a. The large N limit is obtained by assuming that the numbers Nj of 
fermions tend to infinity with equal velocity. At the leading order of the 1/N expansion the RG 
flow is trivial, in the sense that the UV and IR fixed points are indistinguishable. A nontrivial 
RG flow appears at the level of 0(1/N) corrections. In the simplest version of these models, 
the fixed points of these flows coincide with the UV fixed points of the four-fermion models 
(2|, |3[ plus decoupled massless fermions. In a variant of these models there are also fixed points 
of new type. 

The second class of models is a scalar-field counterpart of the first class. Sets of scalar 
fields interact by means of a cr-field. Due to the radiative corrections, ^-potential terms are 
generated. The RG structure of these flows is more involved and the fixed points include known 
conformal field theories (Wilson-Fischer, yj| fixed points) and new ones. 

The theories satisfy certain duality symmetries, associated with the exchange of fermions 
(or scalars) belonging to different sets. The duality symmetries are exact at each order of the 
1/N expansion. Nevertheless, due to some subtleties of the regularization technique ||], the 
implementation of the duality symmetry at the regularized level needs to be studied carefully. 

The integration of the RG flow allows us to write all correlation functions exactly in g to 
a given order in the 1/N expansion. I stress again that the expansion is not around one fixed 
point. Normally, an infinite resummation is necessary to reach an interacting fixed point from 
a Gaussian fixed point. Here, instead, the entire RG flow is expanded at the same time in 1/N 
and no expansion is made in the running couplings g. At each order of the 1/N expansion the 
fixed points and the flow are determined with equal and consistent precision. Typically, both 
the UV and IR fixed points are interacting. 

The large N expansion is crucial to make the integration of the RG flow possible as described. 
To the order 0(1/N), the integration of the RG flow is equivalent to the resummation of the 
leading logs, but in general (in QCD and non-Abelian Yang-Mills theory with matter, for 
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example) the resummation of the leading logs is not sufficient to connect the UV and IR fixed 
points, or the high- and low-energy limits of the theory, although it gives a sort of result beyond 
the perturbative expansion. 

A known theory exhibiting some of the features of the models studied here is the large N 
^93 theory. As a classically conformal (tricritical) theory, it is IR free and its (P(l/iV)-order beta 
function has a nontrivial ultraviolet fixed point Q. The cp® theory, however, does not have a 
duality symmetry. Viewed as a particular case of the models analysed in this paper, it belongs 
to a more general family of theories which do satisfy duality. 

In the study of critical phenomena, stability puts some limits on the values of the ip® coupling 
Q and the UV fixed point of |7[], which lies in the instability region, is superseded by a different 
fixed point. No attempt is made here to apply the approach of ||] to the models of this paper. 
This study is left to a separate research, in view of the possible interest of these theories for 
condensed matter physics. In this paper I study only the tricritical domain. 

Other types of "exact" results about RG flows exist in the literature. For example, in four- 
dimensional supersymmetric theories it is often possible to write an exact beta function The 
NSVZ beta function, however, depends on the subtraction scheme and the present knowledge 
of this scheme does not allow us to integrate the RG flow. The NSVZ formula can be used to 
obtain exact results about the fixed points ||, 

Since the models treated here are classically conformal, it is possible to work with a clas- 
sically conformal subtraction scheme. In practice this means that, by default, the linear and 
quadratic divergences are subtracted away, i.e. the associated renormalized dimensionful cou- 
pling constants are set to zero. This is automatic in the framework of the dimensional regular- 
ization. I recall that in these theories the naive dimensional regularization has to be modified 
adding an evanescent, RG invariant, non-local term to the lagrangian 

The RG flow is scheme independent to the order 0(1/ N). Scheme dependences appear at 
the order 0(1/N 2 ), but scheme-independent correlation functions can be written to arbitrary 
order in the 1/N expansion. 

The models of this paper have been inspired by the study of the irreversibility of the RG 
flow Q. This research singles out, in even dimensions, the specialness of classically conformal 
theories. In odd dimensions the investigation of the irreversibility of the RG flow is more 
difficult, since we miss the tools provided, in even dimensions, by the embedding in external 
gravity. It is useful to construct a large class of classically conformal theories as a laboratory to 
test some ideas on this issue. The interpretation of [10| suggests a more physical choice for the 



expansion parameter, which is the scheme-invariant area of graph of the beta function between 
the fixed points, defined as 

Aa' = J d 3 x|a;| 3 (e(x)e(0)), 

where is the trace of the stress tensor. The study of the quantity Aa' is left to a separate 
work [11], where the implications on irreversibility are discussed. 

In section 2 the fermionic models are studied in detail. In section 3 correlation functions are 
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computed explicitly and the RG interpolation between the fixed points is studied. In section 
4 I solve a particularly simple, but non-trivial, limit, in which the ratios Nj/Nq are sent to 
zero after the large N expansion. In section 5 I study the models with two-component spinors, 
which have peculiar features and fixed points of new type. Section 6 is devoted to the study 
of the scalar models, their beta functions, anomalous dimensions and flows. Section 7 contains 
the conclusions. In the appendix I discuss some subtleties concerning the implementation of 
the duality symmetry in the framework of the modified dimensional-regularization technique of 
ref. 0. 



2 Fermion models 

The model we study was defined in ref. H and describes K + 1 sets of (No, . . . , Nk) fermions 
interacting by means of a <r-field. The lagrangian is 

K Ni 
7=0 i=l 

To simplify various arguments, I first choose four-component complex spinors. The generaliza- 
tion to models with two-component spinors, which have fixed points of new type, is discussed 
later. I work in the Euclidean framework, where ip and ^ are independent from each other. 
The gamma matrices are Hermitian and if) = V^73- A basis with 7^3 = 7^3 and 7J = —72 
is used for charge conjugation (see below). The matrices 74 and 75 of the four-dimensional 
Dirac algebra act on our spinors as flavour matrices. Since the fundamental spinor in three 
dimensions is real (Majorana) and has two components, a set of Nf massless four-component 
complex spinors has a flavour symmetry 0(ANf), while a set of Nf massive four-component 
complex spinors has symmetry 0(2Nf) (g) 0(2Nf). For other details see for example [Hj. The 
field a is a sort of dynamical fermion mass, but our fermions are otherwise massless. The a 
expectation value is zero in the classically conformal phase. 

Symmetries. The theory is manifestly invariant under reflection positivity (Hermitian con- 
jugation combined with a sign-change of the "time" coordinate x%), parity, charge conjugation 
and time-reversal. The P, C and T transformations are 



Pi: 


Xl -> 


-Xl, 


X2,3 


-> £2,3, 


C: 


£1,2,3 


£1,2,3, 


% 


'4 -> 72 


T: 


x 3 -> 


-X3, 


Xl,2 


-» £1,2, 



-> 7174-0/, 



<J — > (7 



a — > a 



i4 7374^/, ipl -> ^[7473, a ^ a. 

Further symmetris are obtained replacing 74 with 75. It is understood that, unless otherwise 
stated, the transformation of ifi is obtained from the one of ip as if were the adjoint of if). In 
the list above, only time reversal requires to transform ip and ^ independently. 
The theory is also invariant under the chiral transformation 

id -> 74^f, <y -o, (2.2) 
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and the analogous transformation with 74 — ► 75 . A consequence of chiral invariance is that the 
graphs with an odd number of external <7-legs vanish. In particular, no vertex <r 3 is generated 
by renormalization. 

The global symmetry group of the theory ( 
the fixed points as discussed below. 

Further simmetries are the g — > 1/g dualities, which I present after fixing some other 
notation. 

If we perform a transformation 

ipf -> 74^, o- -» -a, (2.3) 

for some I only, we change the sign of the corresponding coupling Aj. So, we can assume that 
all of the Xjs are non-negative. 

Renormalization. The field a has no propagator at the classical level, but has a propagator 
proportional to at the quantum level, generated by the fermion bubble. The theory is 

defined using the Parisi large N approach of ref. [|J. The field a has, by definition, the same 
dimension as the fermions, namely (d— l)/2. 

Since the theory is classically conformal, we can choose the classically conformal subtraction 
scheme, in which the linear divergences are just ignored, which is equivalent to say that the 
renormalized fermion mass and the renormalized expectation value of a vanish. 

The dimensional regularization (d = 3 — e) and the minimal subtraction scheme are used. 
Having chosen doublets of complex spinors, we can set the trace of an arbitrary odd number 
of gamma matrices to zero. The symmetries listed in the previous paragraph are manifestly 
preserved at the quantum level. 

It was shown in |IJ that the naive dimensional-regularization technique does not work prop- 
erly, because the fermion bubble generates a a-propagator proportional to 1/ \k 2 ) 2 ~ d l 2 , which 
is responsible for the appearance of T[0]s at the subleading orders. To properly regularize the 
theory, it is necessary to manually convert the a-propagator to l/y/W. This can be achieved 
adding an RG-invariant evanescent non-local term to the lagrangian. 

We choose a preferred set of fermions, say I = 0, as a reference for the large N expansion. 
We define A = Ao, N = Nq, tpi = ip®, etc. We take N large, with rj = Nj/N , aj = Aj/Aq and 

A 2 iV = 4 (2.4) 

fixed. The numerical value in ( p.4| ) is conventional. ([^4|) holds only in d = 3, because for 
d = 3 — e the dimensionality of a and the identity Ab = A// £ / 2 give an evanescent A beta 
function, equal to 

A = 

The couplings aj are exactly dimensionless. We define their renormalization constants by 
(ckj)b = otiZj. The renormalized lagrangian reads 

N K Nj 

c = y: ztfM + \fftzy*a)ik + EE z 0i + A////% 1/2 i?yv)^. 
i=i 1=1 1=1 



2.1) is 0f =o (0(2iV» <g> 0{2N£)), enhanced at 
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Dualities. The duality symmetries of the theory are related to the exchanges of fermions 
belonging to different sets pj. First of all, the theory is evidently symmetric under the exchange 
of 7 with J with 7, J > 0, which means ipj «-> ipj, aj «-> aj and rj «-> rj. 

Let us now exchange the 7th set of fermions, I > 0, with the zeroth set of fermions. The 
theory is symmetric under the duality transformation 



il>T++ifo, ipj^ipj, Nj^Nq, Nj^Nj, A c 



A 



o 



A/ _ 1 aj 1 r ^ K \ 

o^o—^rj, aj — , aj -> — , r j -► — , rj — , (2.5) 
Aq aj aj rj rj 



with J 1,0. The normalization condition (2.4) is duality invariant. The transformation (2.5) 



will be called "the 7th duality transformation" . 

I now prove that, without loss of generality, we can take all of the as to have values between 
zero and one. Suppose that oi\ > 1. (If a\ < 1 then skip to 02-) The first duality transformation 
changes a\ to = l/a% < 1 and ai to a'j = ajja\ for / > 1. Now suppose that a' 2 > 1. (If 
a' 2 < 1 then skip to a^.) The second duality transformation changes a'± < 1 to a'{ = a^/a^ < 1, 
a' 2 to a 2 ' = l/a 2 < 1 and a'j to ct[ = a'j/a 2 , for / > 2. The important point is that the duality 
transformation sends couplings with values < 1 into couplings with values < 1. Clearly, after 
at most K steps we have < ai < 1 for every /. 

In conclusion, it is not restrictive to assume 

< on < 1, V/, 

which will be understood henceforth. 

Fixed points. The conformal fixed points of the RG flow are the models where (aj, . . . oik) 
is any string made of 0s and Is. These theories are always the direct product 

<$ Nf ®FX NlT , (2.6) 

of a theory ^N f of Nf free massless fermions and the interacting conformal field theory (called 
on a in ref. 0) identified by the lagrangian 

N„ __ 

C = J2^ + X °)^- ( 2 - 7 ) 

For practical convenience, this theory will be called FT,^ henceforth. Here Nf is the total 
number of decoupled fermions (aj = 0) and N a is the total number of fermions coupled to a with 
ai = 1. The theory (2.6) has an enhanced flavor symmetry group 0(4JVf) <8> 0(2N a ) (%> 0(2N a ). 



The conformal field theory FT, jv ct coincides with the UV fixed point of the four-fermion models 



N - 1 ( N - V 

i=l \i=l / 
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(a) (b) 




(C) (d) 

Figure 1: Order 0(1/N) graphs 

Following Q these conformal field theories can be formulated per se, in a classically conformal 
framework, with no reference to a flow originating them as its large- or small-distance limit. 

The zeroth order of the large N expansion. The cr-propagator is dynamically gener- 
ated by the one-loop fermion bubble (graph (a) of Fig. 1) and reads 

1 K 

(a(k)a(-k)) = —D(a J: rj), D' x {aj,rj) = 1 + £ ajrj. (2.8) 



To the leading order of the large N expansion the quantum action reads 

K Nj , A , 

r = J d d *££^#f + J j^\*(k)\ 2 VtfD-\aj,rj), (2.9) 



7=0 i=l 



This is a Gaussian theory, with a non-local kinetic term for a. The couplings are inert and 
no flow takes place. In this limit, the UV and IR fixed points are indistinguishable, and the 
constant D(aj,rj) can be reabsorbed in a. It is immediate to check that the zeroth-order 
expression (|2.9|) of T is invariant under the duality transformation (|2.5[). 



The first order of the large N expansion. The RG flow can be inspected computing the 
0(1/N) corrections. The relevant graphs are the one-loop fermion self-energy and the one-loop 
vertex (see Fig. 1) nm. A two-loop graph (d) containing a fermion triangle could in principle 
contribute to the vertex, but it is zero because of the conservation of chirality. 

After resumming the fermion bubbles, the number of relevant graphs becomes finite at each 
order of the 1/N expansion. This means that the a/-dependence can be determined exactly, 
order by order in the 1/N expansion. Moreover, this knowledge is sufficient to connect the fixed 
points, so that it is possible to write a/-exact beta functions and correlation functions to any 
given order of the 1/N expansion. The order-0(l/iV) expressions of the beta functions and 
anomalous dimensions allow us to resum the leading logs (see next section). 
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The order-0(l/iV) renormalization constants can be written using the results of [fl| and 



read 



J V>' 



1 



-aj/i £ D(aj,rj) 



3tt 2 Ne 
16 

Zi = l + 37T 2 Ne ( a/ ~ 1 V~ £ £>(aj, »v) 



1 + 



16 



3vr 2 iVe 



jix e D(aj,rj), 



(2.10) 



These expressions are correct also for 1 = 0. The factors D(aj,rj) are brought by the a- 
propagators. The factor fi~ £ should be understood as a short form for AgiV/i -e /4, so that, 
using (HD, 

1 / AjiV^ y 



lim ^ 
e^o dln/i 



-to. 



We immediately derive anomalous dimensions and beta functions: 
ldlnZj, 1 , , ldlnZ^ 



Pilot) 



2 dln/i 3vr 2 iV 



a I D(aj,rj) 



la 



2 dln/i 



3vr 2 iV 



D(aj,rj), 



16 



dlnZj 

din// ~ 3^2iV 



ai(ai - l)D(aj,rj). 



(2.11) 



The duality (2.5) can be checked explicitly in these formulas. The only non-straightforward 
step in this check is the trasformation of j a : 



la -> la 



1 

2o~ 



-Pi 



ll 



(2.12) 



which is implied by a — ► a^Jrj\j/\o. We are going to check duality also in the exact expressions 



of the correlation functions of the next section. From the duality relation ( 2.12 ) we recover the 
formula for the beta function derived in Ififl, 



which is exact. 



Solution of the RG equations. The RG flow can be integrated immediately in the space 
of couplings. We have 

daj daj 



ai(ai - 1) ajiaj - 1) 
for every I and J. The solution of ( |2.13| ) is 

i--i = c„(-L-i 



(2.13) 



(2.14) 



with no summation over J. The matrix Cjj is a matrix of arbitrary constants. Restricting to 
< oli < 1, self-consistency of ( |2.14 ) demands that Cjj satisfies 

1 



Cu > 0, CjjCjl — Cil, 



c 



I J 
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where, again, no summation on J is understood. The most general solution to these constraints 
is determined by the constants Ci = Cu with Cjj = Cj/Cj. We can write 

«/ = ..ft* 1W (2-15) 
1 + (pi — l)a 

where a = a\. The choices Cj = 0, 1, oo for some Is correspond to freeze some sets of fermions 
or glue them to other sets of fermions. 

The dependence on the energy is obtained solving the RG equation 

16 a(a — 1) da 

(0) = 3^7V 1 + ri a + E?= 2 Cjrja/ (1 + (Cj - l)a) = ~ dt ' 

where i = ln|x|/u and \x\ is the overall scale of the correlation function. This equation can 
always be solved implicitly in the form t = t(a). In general, the function t(a) can be inverted 
only numerically, but in the simplest case (K = 1, r\ = 1) it can be inverted also explicitly. 
The solution is then 

«(VM) = i + u - v / ^7^), u = (1 (i^r 16 /^) . 

We have a(l/|a;|) — * for |x| — > (ultraviolet) and a(l/|a;|) — > 1 for |x| — > oo (infrared). For 
K = 1 and ri = r generic we have 

(|duT 16/(37rSjv) - f 1 ~ "^/'^^ " " ' 



a(l/|x|) V 1 — "0-0 

3 Correlation functions in fermion models 

In this section I show how the integration of the RG flow allows us to write all correlation 
functions exactly in a to any given order in 1/N. Here it is understood that ln/i ~ iV, because 
the powers of (ln/x)/iV are resummed. In particular, to write the correlation functions at the 
leading-log level, it is sufficient to know the beta function and renormalization constants to the 
order 1/N. For concreteness, I choose K = 1 and r = 1. This model has two sets of fermions, 
ijp and ip , and one coupling constant a. 

The a two-point function. Using the Callan-Symanzik equations we know that the a 
two-point function has the form 

<a(*M0)> = Z-\a{l/\x\)M^)) A{ ^J^ ) ■ 

At the leading-log level, we need 70- to (D(l/N), given by (2.11), and A to 0(1), which is just 
1/(1 + a), from Q. We find 

Z a (a{l/\x\),a{n)) = exp -2 / 7<j(aO d 111 ^ " 



ln/i ' ' / "(l/l^l) 1 - 
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and therefore 

(a(x)a(O)) 



1 a(l/\x\) l-a(ji) 



1 



1 1 — a{jj) 



(3.1) 



2vr 2 |x| 2 l-a 2 (l/|x|) a{ji) An 2 \x\ 2 ^u{u + 2) a(p) 
It is amusing to check the invariance of this correlation function under the duality transforma- 



tion ( |2.5D . In doing this, remember that two as get a factor a(fi) in the transformation. 

The exact formula ( |3.1[ ) exhibits the dependence on the reference scale /j,, which survives at 
both critical points and is different in the UV and IR limits. The anomalous dimension varies 
from 1 — 8/(37r 2 iV) in the ultraviolet, to 1 — 4/(37r 2 iV) in the infrared, as expected. The UV 
and IR behaviors are 



respectively. 



1 



16/(3tt 2 A0 



27r 2| x |2(l-8/(3^7V)) I-q^) 



and 



8/(37r 2 AT) 



4^2^12(1-4/(3^^)) 



The fermion two-point functions. With the same procedure, the fermion two-point 
functions are found to be 



<?V)v> u (o)> 



(x)^(0)) 



i 



47T | X | 3 

i 

47r|x| 3 



1 — a(l/|x|) Oi(fi) 

a(l/\x\) 1 — a{(j) 
l-a(l/\x\y s 



1 — a(/x) 

There formulas are interchanged by the duality transformation (|2.5| ) (observe that under duality 
a(l/|x|) — > l/a(l/\x\) and a(/x) — > l/a(/j,)). 

The cr 2 two-point function. As a further example, we study the two-point function of the 
composite operator a 2 . This operator has dimension 2, but to the order 1/N it does not mix with 
the fermion bilinears ip if) and ip ip 1 (it is sufficient to consider only one such bilinear, since the 
a field equation relates the two to each other). The calculation of the order-0(l/iV) anomalous 
dimension of a 2 requires the evaluation of the divergent parts of two two-loop diagrams with 
a fermion loop (see Fig. 2). Since both the fermions with / = 0, 1 circulate in the loops, the 
result is multiplied by 1 + a 2 . Since there are two cr-propagators in these diagrams, the result 
has also a factor 1/(1 + a) 2 . We obtain, taking into account of the combinatorial factor and 
the number of ways in which the fermions can circulate in the loop, 



and therefore 



Z„2 Zo 



1 + 



Nir 2 e 



-2c- 



1 



1 + a 2 
(1 + a) 2 

2a + 3a 2 



3Nir 2 (1 + a) 2 ' 
Integrating the RG equations we find 

A 2 (a(l/\x\)) ^u{u + 2) 



(a 2 (x)a 2 (0)) = Z- 2 2 (a(l/\x\), «(/,)) 



a{n) 



4vr 4 |x| 4 (l-a 2 (/u))(l + a(^)) 



2 ' 
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Figure 2: Renormalization of the operator a 2 



The function A 2 {a) can be read from the order-O(l) expression of the correlation function, 
which is the a-bubble. We have A2(a) = 1/(1 + a) 2 . 

Again, the correlation function transforms correctly under the duality (2.5), 

In the ultraviolet and infrared limits we have 

1 /i- 16 /( 3 ^)(l-a(/i)) , 1 /i- 8 /( 3 - 2iV )v^(M) 

and 



87r 4| x |2(2+8/(3^)) (l+a(^)) 3 4vr 4| x |2(2+4/(3^)) (1 + a (^)) 3 

respectively. 

In more generality, the anomalous dimension of the operator a 2 is 

K 

* ,-,■>. , / , 

7<r2 



:D 2 (aj,rj) (l + J2ajrj(3aj - 2)^ 



3ir 2 N 

and its transformation under the Jth duality reads 7^2 — > 7^2 — (3j/aj. 

Other correlation functions. The integration of the RG flow solves the a-dependence 
and allows us to write the correlation functions exactly in a to any given order in 1/N. In this 
expansion, if s denotes the overall scale of the correlation function, lns/i is considered of order 
N and the powers of (In sjj)/N are resummed. In the previous paragraph I have considered 
two-point functions explicitly, but the procedure is general. 

Let us consider the three-point function 

(a 2 (x)a 2 (y)a 2 (0)) = G(\x\,\y\,\x-y\). 

This correlation function can be more conveniently expressed in terms of the variables 

Ixl \y\ (x — y) 2 

s = \x + \y + \x - y , u 1 = — + —, v 2 = ————. 

\y\ \x\ \ x \\y\ 

The variables are scale-independent and s is the overall scale. We can write 
G{s,v 1 ,v 2 ,a(n),N) = — ^ N k fk{vi,v 2 ,a.(n),N). 

k 

For the leading-log approximation it is sufficient to know the functions to the leading order. 
The Callan-Symanzik equations give 

1 _ 3 

G(s,ui,U2,a(n),N) = ~^Z a 2 {a(l/s),a{fjL))G Ted (u 1 ,U2,a(l/s),N) 
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and it is sufficient to compute Gred^ii a(l/s), N) to the leading order in 1/N. We have 

(^Wivwm - ^ |x|b ; |x _ y| w.)w.) + 2» 3 " ( (i _ a 2( A i°)(i +a ( M ))2 ) 3/2 ' 

where u(s) = (s/j)~ 16 /( 37r2Ar )(l — a(/j)) 2 /(2a(/i)). If we choose a different definition of the 
variable s, formula (|3.2j) is unaffected at the level of leading logs. 



4 The r = limit 

A particularly simple limit is the case in which the rjs are zero. Setting rj = 0, I = 1, . . . K , 
after expanding in 1/N is not equivalent to kill the fermions with / = 1, . . . K. The large N 
and 77 — > limits conflict and leave a non-trivial remnant. When 77 = 0, however, duality is 
not manifest. We can see this immediately from the expression of the beta function: 

16 

The RG flow interpolates between a = and a = 1, but cannot reach the fixed point at a = 00. 
To reach a = 00 we have to perform a dual limit in which some or all of the 77s tend to infinity. 
The solution reads 

"/(VM) - 



a(fj,) + (1 - a(fi))Cy\\x\fi)- 16 /^ 2N ) ' 
where a = a\ and C\ = 1. Formula fl2.15|) still holds. The anomalous dimensions are 

j 1 8 8 

^ ~ 3^N ai ' 7<T ~ ~3^iV' 7,72 ~ 3^iV' 

7o- and 7^2 do not depend on the energy. The other results can be generalized straightforwardly 
to this case. 

5 Models with two-component spinors and new fixed points 



Let us now consider a variant of the model (|2.l| ) in which the spinors have two complex com- 
ponents instead of four. 

We have no chiral invariance ( |2.2[ ) and the transformation ( |2.3|) does not apply. For this 
reason, the signs of the couplings gj cannot be assumed to be positive. The P and T transfor- 
mations are modified as follows: 

Pi: xt^-xt, x 2 ,3^^2,3, ipi 7i^f > cr -> -a ; (5.1) 
T : x 3 -> -x 3 , xi >2 -> xi,2, V'j 7 73^, ^73, cr -o - . 



The theory is still invariant under the duality (2.5) and we can assume that —1 < gi < 1. 
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To avoid the problem of extending to d dimensions a Dirac algebra in which the trace of an 
odd number of gamma matrices does not vanish, it is better to keep d = 3 and regularize the 
theory modifying the ip- and cr-quadratic terms with a cut-off as follows: 

H (l " |f) V>, and l\a(k)\ 2 Vtf (l + ^ D(gj,n). 

This regularization framework manifestly preserves the P and T symmetries Q5.1| ) (a Pauli- 
Villars framework does not) and ensures that no a 3 term is generated by renormalization. The 
graphs with an odd number of a external legs are finite, but need not vanish. 

The wave-function renormalization constants and beta functions are the same as in ( |2,10D 
and ( |2.1l| ) with iV — > N/2. In particular, the gi beta functions read 

^-f^^^-W-o). (5.2) 

We see that the fixed points of this model are {gi}= any string made of 1,0, —1. These fixed 
points have the form 

where Nj is the number of decoupled fermions (coupling gj = 0), N a is the number of fermions 
interacting with a with coupling gj = 1 and N' a is the number of fermions interacting with a 
with coupling gj = — 1. The conformal field theories i^^A^A^ are defined by the lagrangian 

i=l i=j 

and have symmetry 0(2N (7 )^)0(2N' (7 ). The subscript 2 in i^Sjv^jv' means that the spinors are 
two-component. In general (r 7^ 1), the conformality of ( |5.3| ) does not seem to be implied by 
more general principles than the vanishing of the beta function (|5.2| ) at g = — 1. The symmetries 
of the model allow us to construct a local, dimension-3 operator O that does not vanish using 
the field equations (for example aipijj) and could appear in the trace anomaly, = j30. Instead, 
at r = 1 the two-component spinors ip and x can be grouped into a four-component spinor and 
we have i^Sjv.jv = ^kSjv, whose conformality follows from the general arguments of [||. 
The i^SjVo-,^/ theories have also a duality 

tp ^ X, a^-a^/r, A — > A/a/t, r -> l/r, N a <-> iV£, 

where r = iV;/iV ff and A 2 ^ = 8. 

The order 0(1/ N) anomalous dimensions 7! and 7^ of the theories i^EAT.A" coincide with 
those of the theories i^^AM-A" (and ^4^^/2,^725 F^(N+N')/2)t but the symmetry groups differ. 

It is straightforward to extend the other results of the previous sections to the flows inter- 
polating between the F~En,n' models. 
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6 Scalar models 



1 K Nj 



We now consider the scalar models 

K Nj 

z /=0 i=l 

with symmetry group <S>f=o O(Nj), possibly enhanced at the fixed points. The A/s are not nec- 
essarily positive and V%{ip) denotes a potential term of the form <£> 6 induced by renormalization. 

Scalar Sjy conformal field theories. I define the scalar Y,jy conformal theories by means 
of the lagrangian 

1 N 

C = -Y J [{d^? + ih<^) 2 ]: (6.2) 

i=l 

which I briefly denote with STi^. This theory coincides with the UV (Wilson-Fischer) fixed 
point of the O(N) sigma model. In the approach of Q, the SY,n model can be studied per 
se, with no reference to an RG flow generating it as a critical limit. The construction of the 
theories n follows the lines of the construction of the fermionic theories FT, n of JTJ . A caveat 
is that renormalization generates a (p e term. Due to the O(N) symmetry, this term can only 
have the form 




(6.3) 

=i / 

and so can be reabsorbed with an (imaginary) translation of the field a. This is equivalent to 
say that ( |3.3| ) is "trivial" because proportional to the a field equation. The a field equation 
formally sets the composite operator if 2 to zero, but does not trivialize the theory. More details 
on the use of the a field equation can be found in the appendix. 

The theory is classically conformal and we choose the classically conformal subtraction 
scheme, in which the renormalized (p 2 - and y? 4 -terms are set to zero by default. 

The a-propagator is generated by the one-loop scalar bubble, which is proportional to 
{k 2 ) d / 2 - 2 . If unmodified, a er-propagator ~ (k 2 ) 2 d ^ 2 generates T[0]s at the subleading orders. 
It is therefore necessary to add an evanescent, RG-invariant, non-local term to the lagrangian 
and convert the a-propagator to \/W. This means, in particular, that the field a has dimension 
(d + 1) /2 and Ab = A/// 2 . The 0(1/N) diagrams are the same as in Fig. 1. The diagram (d) 
does not vanish identically, but is finite. Finally, the 0(l/iV) results for and 7 CT are equal to 
those of the FY,n models ( |2.7| ), with the replacement N — > N/4. 



The running models. Let us now consider the theories fl6.ip . We can still use the 

classically conformal substraction scheme, and therefore assume that the terms of the form if 2 
and <£> 4 are subtracted away by default. However, nothing forbids the generation of terms of 
the form p e and this time they cannot be removed with a translation of the field a. The reason 
is that a multiplies 

f>(/) 2 , 

7=0 
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but the ip divergences can be proportional to 

(f>?V) 2 ) (l>IK^) 2 ) (e^) 2 ) 

with m, n and p generic. It is therefore necessary to introduce the most general y? 6 -vertices in 
the lagrangian, multiplied by independent couplings. 

For simplicity, I focus on a model with two sets of scalar fields, tpi, i = 1,...N and Xji 
j = 1, . . . Nr, K = 2. The lagrangian reads 

(6.4) 

The factors in front of the Vg-vertices are chosen for notational convenience. 
The duality transformation is 

g -, ip <-> x, r->-, A 60 A 6, A -p, cr -> cr^y^- (6.5) 

I have dropped vertices of the form ip 4 x 2 and V^X 4 ; since they can always be reabsorbed 
with a cx-translation while keeping duality manifest. The remaining ^-vertices are redundant 
(one suffices), but cannot be reduced to a single vertex in the full space of parameters using 
the ex-field equation. We have to divide the parameter space into two charts, one contaning 
g = and the other containing g = oo. In the first chart, the if 6 vertex is trivial and can 
be reabsorbed into a er-translation. However, the x 6 vertex cannot be reabsorbed into a a 
translation at the point g = 0. Therefore, in the first chart we remove the ip G vertex setting 
Aeo = and take Ao6 as the unique ^-coupling. In the second chart, which contains g = oo, we 
have a symmetric situation: the x 6 vertex can be reabsorbed away, but the ip 6 vertex cannot, 
so we set Ao6 = and the true coupling is Ago- 

In each chart we have an appropriate specialization of the duality transformation (|5.5|). For 
example, when we apply duality in the first chart, we transform the potential \§Qtp & /N 2 into 
A6oX 6 /(N 2 r 2 ). Then, we can use the cr-field equation x 2 = — V 2 /g and rewrite the potential as 
— Xsotp 6 / (r 2 g 3 N 2 ) . To recover the initial potential we define Xqo = —r 2 g 3 Xeo, so that duality in 
the first chart reads 

g -, ¥> X, r -, A 60 -» -rVAeo, A -» — , a -> ag^fr. (6.6) 
g r yr 

A symmetric argument allows us to specialize the duality transformation to the second chart. 
Invariance under duality is a powerful tool to check the calculations. 

The leading and subleading order in the cr-sector. The effective cr-propagator is 
obtained resumming the scalar bubbles and reads 

(a(k)a(-k)) =VtfD(aj,rj). 
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To the order 0(1/N), the </? 6 and x 6 vertices have no effect on the renormalization of the 
scalar kinetic terms and the cr-scalar-scalar vertices. This can be proved observing that the 
relevant graphs are in one-to-one correspondence with the graphs of Fig. 4 contaning four 
scalar external legs exiting from the same six-leg vertex. I call these graphs Q^2 for clarity. 

i) The graphs contributing to the renormalization of the scalar propagator are obtained 
from C/4 2 suppressing the four external scalar legs exiting from the same six-leg vertex. We 
obtain the graph (b) of Fig.l plus tadpoles, finite graphs and a graph with a scalar bubble (from 
the third of Fig. 4). This graph should not be counted, since the scalar bubbles are resummed 
into the a-propagator. 

ii) The graphs contributing to the renormalization of the cr-scalar-scalar vertices are ob- 
tained from replacing the four external scalar legs exiting from the same six-leg vertex 
with one external a leg. We obtain the graphs (c) and (d) of Fig. 1 plus tadpoles and finite 
graphs. 

In summary, the order 0(1/N) anomalous dimensions and beta functions are exactly the 
same as m fl2.1l| ), with N -> N/4. In particular, we have 

32 g(g 2 -l) 

Pg 



3Ntt 2 l + rg 2 
The zeroes are at g = — oo, —1, 0, 1, oo. 

The (y9 6 -terms. Now we study the induction of a V6-coupling by renormalization. The 
order 0(1/ N) beta function of the y? 6 -coupling reads in the first chart 

Pm = 3^2(1 + ^2)3 [-(! " 9f + 3A 60 (l + 3g 2 + 2rg 2 + 6r 5 3 + 4rV 

-27A| r 5 3 (l - r V) - 27A 3 rV] . (6.7) 

Graphs and details on the calculation are given in the appendix. This formula transforms 
correctly under the duality (|6.6| ), i.e. 



Pea -» -rV/%o - 3A 60 r 2 5 2 /? ff . (6. 



In the second chart we have 
32 



A)6 = 3A r r7r 2 (1 + rff2) 3 [ r V(l - 9? + 3A 06 r 5 2 (4 + 6rg + 2rg 2 + 3r 2 g 2 + r 2 g l 



+27A 2 6 (l-rV)"27A 3 
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Fixed points. At g = oo and g = we work in the first and second chart, respectively, 
and find the beta function of ref.s Q, namely 

288 288 
/%o = ]y^2^6o(l ~ ^6o)> A)6 = ^y^ 'W 1 ~ ^oe)- (6-9) 
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The solutions of /3go = are Ag = and Ag = 1. At Ag = the fixed point is the direct product 
of an 5Sjy r model and N free scalars. At Ago = 1 the fixed point is the direct product of SHff r 
and the interacting ip® fixed point. In both cases, the flavor symmetry group is O(N) ® O(Nr). 
The negative values of Ago and Ago are unphysical at g = oo and g = 0, becuse the potential is 
unbounded from below. 
For g = 1 we have 

= Wa+r)3 I 4(1 + " 9A60r(1 " " 9X ^ 3 
The solutions of /3go = are Ag = and the two interacting fixed points 

1 + r 



A 60 



6r 2 



(3(r - 1) ± y/9-2r + 9r 2 ) , 



which are interchanged by duality. At A60 = the fixed point is an 5£jv(i+r) model. 
For g = — 1 we have 

32 r 

&o = 3iV7r 2 (1 + r) 3 [" 8 + 12A 6o(! " r + r 2 ) + 21\ 2 m r{\ + r 2 ) - 27AgV 

The equation /3go = has three real solutions, two positive and one negative. For r — > oo the 
largest positive solution tends to 1 from above, the second positive solution Ago tends to + 
and the negative solution tends to CP. 

The fixed point (g, Ago) = (0, 1) is UV unstable, while the fixed points (1,0) and (— 1, Ago) 
are IR stable. The other fixed points are partially stable and partially unstable. 

The RG flows are plotted in Fig. 3. Since at g = the negative values of Ago are unphysical, 
we can argue that the unphysical region is the region below the line connecting the three fixed 
points with Ag < 0. 

A particularly simple case is the limit r — ► oo, where /3go is independent of g and coincides 
with (|6.9|). Here the RG flows of g and Ago separate. 



7 Conclusions 

In the realm of high-energy physics, classically conformal theories play a special role. Taking 
inspiration from massless QCD, where hadrons and glueballs acquire masses dynamically, it is 
tempting to think that the ultimate theory of the universe contains no dimensionful parameter 
of classical origin, the classically conformal subtraction scheme should be universally used and 
the dynamical generation of a scale /x by dimensional transmutation originates all scales and 
dimensionful parameters of nature, including the Newton constant. 

Another context in which the specialness of classically conformal theories is singled out, at 
least in even dimensions, is the irreversibility of the RG flow. It would be interesting to know 
if the RG flow is irreversible also in odd dimensions. The flows of this paper are used for this 
investigation in ref. |11]. 
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Figure 3: The fixed points and RG flows of the scalar models 



It is useful to have a large class of classically conformal models to test ideas around the 
issues just mentioned and the strongly-coupled limit of quantum field theory. I believe that 
the models studied in this paper can help us answering various unsolved problems and have a 
number of potential applications, maybe also in condensed-matter physics. 

I have explored two classes of theories, describing sets of fermions and scalar fields interacting 
by means of a d-field, with a variety of fixed points, known and new. Typically, both end 
points of the flows are interacting conformal field theories. The remarkable properties of these 
theories are that i) at each order of the 1/N expansion the RG flows can be integrated for 
arbitrary values of the running couplings; ii) "exact" beta functions, anomalous dimensions and 
correlation functions interpolating between the UV and IR conformal limits can be written; in) 
the models have manifest g — ► 1/g dualities at each order of the 1/N expansion. Duality puts 
severe constraints on the formulas and is a powerful tool to check the calculations. In special 
limits (r — > in the fermionic models, r — > oo in the scalar models) the integration of the flow 
equations further simplifies, yet the flow remains non-trivial. 

Appendix: diagrammatics and duality in the scalar models 

The graphs contributing to the renormalization of (p G to the order 0(1/N) are shown in 
Fig. 4. Some diagrams have obvious counterparts with two, four and six external x~^ e E s - 
The divergent terms proportional to (p 2 x 4 i V 4 X 2 an d X 6 are converted to ip 6 using the a-field 
equation (p 2 + g\ 2 = 0. (For subtleties in the use of this field equation see below.) Several 
diagrams are finite, but have been drawn for completeness. 
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Figure 4: Renormalization of c/? 6 in the scalar models 



The beta function (qJ) is invariant under the duality symmetry ( |6.6|) and (|3^). However, 
the renormalization constant Zgo is not. The complete expression of XqoZqq is very long and 
contains terms with the first three powers of Aeo plus a term independent of Ago- Here I focus 
on the terms proportional to Aq , for simplicity, which are 



288 ^60 
iVvr 2 



-Is 



3/i 



-Be 



4c 



+ 



-6s 



-7 s 



(7.10) 



5e(l + rg 2 ) 2e(l + rg 2 ) 2 7e{l + rg 2 f)' 

Expanding ^~ m£ /(me) as l/(me) — In// we see that the terms proportional to ln/i, which 
contribute to the beta function, transform correctly under duality, while the pole terms, which 
have no physical significance, do not transform correctly. 

These facts are explained as follows. The symmetry of the theory under the duality trans- 
formation ( |6.6| ) in the first chart understands the use the a field equation. Naively, the a field 
equation is ip 2 + gx 2 = 0, but the modified regularization technique of M demands to add the 



non local term 



1 + rg 2 



d d k 



K*0l s 



i 



i 



X 2 B N 



-s/2" 



2 J (2vr) dl " v_/l \ 16 

to the action S. The complete a field equation reads, in momentum space, 



ip 2 (k)+g X 2 (k) 



2i(l + rg 2 ) 1 



A 



1 



A|iV 
16 



(*>) 



-s/2" 



a(k). 



The evanescent right-hand side of this formula provides extra vertices which do contribute to 
the pole terms. However, these corrections are always proportional to 



-mis 



-m 2 s 



VTL\£ 



m 2 e 



(7.11) 
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for some m\ and m,2 and so do not affect the In /i, terms. This fact can be easily verified 
in a number of examples. Moreover, the terms proportional to the a field equation 5S/5a 
cannot be neglected. They give other contributions of the form ( |7. 1 1| ) , which can be easily 
computed functionally integrating by parts. Putting these ingredients together, we discover 
that the mismatch between Zqq and its dual is due to the use of the a field equation and the 
difference between the naive and regularized a field equations. It is not a lengthy work to 
make a quantitative check for the terms ( 7. ICQ . The important point is that the mismatch has 



no physical significance. We can neglect these nuisances if we limit ourselves to check duality 
within the physical quantities. 
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